## Loading required package: foreach
##
## Attaching package: 'foreach'
## The following objects are masked from 'package:purrr':
##
## accumulate, when
## Loading required package: iterators
##
## Attaching package: 'data.table'
## The following objects are masked from 'package:lubridate':
##
## hour, isoweek, mday, minute, month, quarter, second, wday, week,
## yday, year
## The following objects are masked from 'package:dplyr':
##
## between, first, last
## The following object is masked from 'package:purrr':
##
## transpose
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.088587 -0.006528 -0.001660 0.005455 0.060812
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0102102 0.0002809 36.35 <2e-16 ***
## test_set$FCH4 0.1937750 0.0096910 20.00 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01042 on 1800 degrees of freedom
## Multiple R-squared: 0.1818, Adjusted R-squared: 0.1813
## F-statistic: 399.8 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.214399 -0.004156 -0.000870 0.003302 0.066976
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063025 0.0001128 55.9 <2e-16 ***
## result$FCH4 0.5080646 0.0041981 121.0 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008289 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.6703, Adjusted R-squared: 0.6703
## F-statistic: 1.465e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.056812 -0.006415 -0.001517 0.005257 0.045742
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0098642 0.0002803 35.19 <2e-16 ***
## test_set$FCH4 0.2669940 0.0108776 24.55 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01016 on 1800 degrees of freedom
## Multiple R-squared: 0.2508, Adjusted R-squared: 0.2504
## F-statistic: 602.5 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.140559 -0.003996 -0.000804 0.003153 0.091726
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0060996 0.0001033 59.07 <2e-16 ***
## result$FCH4 0.5465918 0.0038445 142.18 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007591 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7373, Adjusted R-squared: 0.7372
## F-statistic: 2.021e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.046781 -0.006555 -0.001258 0.005046 0.062954
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0100890 0.0002652 38.05 <2e-16 ***
## test_set$FCH4 0.2342332 0.0100435 23.32 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009778 on 1800 degrees of freedom
## Multiple R-squared: 0.2321, Adjusted R-squared: 0.2316
## F-statistic: 543.9 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.084343 -0.004111 -0.000991 0.003176 0.124377
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0062639 0.0001049 59.72 <2e-16 ***
## result$FCH4 0.5355601 0.0039053 137.14 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007711 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.723, Adjusted R-squared: 0.723
## F-statistic: 1.881e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.101141 -0.006816 -0.001937 0.005463 0.065007
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0105040 0.0002833 37.08 <2e-16 ***
## test_set$FCH4 0.2200137 0.0101462 21.68 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01049 on 1800 degrees of freedom
## Multiple R-squared: 0.2071, Adjusted R-squared: 0.2067
## F-statistic: 470.2 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.222270 -0.004195 -0.000961 0.003341 0.114902
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063836 0.0001097 58.21 <2e-16 ***
## result$FCH4 0.5234734 0.0040831 128.21 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008062 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.6953, Adjusted R-squared: 0.6952
## F-statistic: 1.644e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.071778 -0.006547 -0.001534 0.005158 0.072791
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0091863 0.0002951 31.13 <2e-16 ***
## test_set$FCH4 0.3050215 0.0112129 27.20 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01072 on 1800 degrees of freedom
## Multiple R-squared: 0.2913, Adjusted R-squared: 0.2909
## F-statistic: 740 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.136430 -0.003859 -0.000902 0.003060 0.069183
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.005580 0.000107 52.16 <2e-16 ***
## result$FCH4 0.580003 0.003983 145.61 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007865 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7464, Adjusted R-squared: 0.7464
## F-statistic: 2.12e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.104647 -0.006499 -0.001978 0.004736 0.054491
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0104210 0.0002751 37.88 <2e-16 ***
## test_set$FCH4 0.2350185 0.0104403 22.51 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.0103 on 1800 degrees of freedom
## Multiple R-squared: 0.2197, Adjusted R-squared: 0.2192
## F-statistic: 506.7 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.223639 -0.004109 -0.001156 0.003031 0.103082
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0064960 0.0001072 60.61 <2e-16 ***
## result$FCH4 0.5328294 0.0039903 133.53 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007879 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7122, Adjusted R-squared: 0.7122
## F-statistic: 1.783e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.101082 -0.006468 -0.001583 0.004925 0.052405
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0104478 0.0002661 39.26 <2e-16 ***
## test_set$FCH4 0.2210462 0.0099385 22.24 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009968 on 1800 degrees of freedom
## Multiple R-squared: 0.2156, Adjusted R-squared: 0.2151
## F-statistic: 494.7 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.224632 -0.004161 -0.001050 0.003044 0.115431
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.006478 0.000107 60.52 <2e-16 ***
## result$FCH4 0.530006 0.003985 132.99 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007869 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7106, Adjusted R-squared: 0.7105
## F-statistic: 1.769e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.065688 -0.006424 -0.001524 0.004955 0.054633
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0101566 0.0002679 37.91 <2e-16 ***
## test_set$FCH4 0.2342388 0.0102115 22.94 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009869 on 1800 degrees of freedom
## Multiple R-squared: 0.2262, Adjusted R-squared: 0.2258
## F-statistic: 526.2 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.165873 -0.004055 -0.001017 0.003211 0.104060
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0062579 0.0001052 59.51 <2e-16 ***
## result$FCH4 0.5378033 0.0039151 137.37 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007731 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7237, Adjusted R-squared: 0.7237
## F-statistic: 1.887e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.052091 -0.006942 -0.001918 0.004887 0.075204
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.009605 0.000320 30.02 <2e-16 ***
## test_set$FCH4 0.305081 0.012653 24.11 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.0115 on 1800 degrees of freedom
## Multiple R-squared: 0.2441, Adjusted R-squared: 0.2437
## F-statistic: 581.4 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.090522 -0.003988 -0.001092 0.002956 0.097211
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0055438 0.0001109 49.98 <2e-16 ***
## result$FCH4 0.5943966 0.0041298 143.93 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008155 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.742, Adjusted R-squared: 0.7419
## F-statistic: 2.072e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.094260 -0.006977 -0.001810 0.005738 0.063195
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0107564 0.0002876 37.40 <2e-16 ***
## test_set$FCH4 0.1842998 0.0096624 19.07 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01089 on 1800 degrees of freedom
## Multiple R-squared: 0.1681, Adjusted R-squared: 0.1677
## F-statistic: 363.8 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.231763 -0.004264 -0.001048 0.003292 0.088957
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063058 0.0001207 52.24 <2e-16 ***
## result$FCH4 0.5282328 0.0044939 117.55 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008873 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.6573, Adjusted R-squared: 0.6572
## F-statistic: 1.382e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.099896 -0.006180 -0.001648 0.005092 0.063468
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0105304 0.0002713 38.81 <2e-16 ***
## test_set$FCH4 0.2076361 0.0098542 21.07 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01009 on 1800 degrees of freedom
## Multiple R-squared: 0.1979, Adjusted R-squared: 0.1974
## F-statistic: 444 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.224931 -0.004125 -0.000962 0.003276 0.085707
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0064549 0.0001093 59.08 <2e-16 ***
## result$FCH4 0.5204512 0.0040680 127.94 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008033 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.6944, Adjusted R-squared: 0.6943
## F-statistic: 1.637e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.066378 -0.006345 -0.001483 0.004949 0.061575
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0102235 0.0002733 37.41 <2e-16 ***
## test_set$FCH4 0.2338459 0.0101743 22.98 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009994 on 1800 degrees of freedom
## Multiple R-squared: 0.2269, Adjusted R-squared: 0.2265
## F-statistic: 528.3 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.164817 -0.003991 -0.000915 0.003311 0.122674
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0062297 0.0001059 58.81 <2e-16 ***
## result$FCH4 0.5325967 0.0039436 135.05 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007787 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7169, Adjusted R-squared: 0.7168
## F-statistic: 1.824e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.081565 -0.006991 -0.001949 0.005499 0.096607
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.010179 0.000324 31.41 <2e-16 ***
## test_set$FCH4 0.278310 0.012328 22.58 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01172 on 1800 degrees of freedom
## Multiple R-squared: 0.2207, Adjusted R-squared: 0.2202
## F-statistic: 509.7 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.174728 -0.004079 -0.001184 0.003038 0.138382
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.005721 0.000114 50.16 <2e-16 ***
## result$FCH4 0.580983 0.004246 136.83 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008384 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7221, Adjusted R-squared: 0.7221
## F-statistic: 1.872e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.045889 -0.006677 -0.001710 0.004986 0.066562
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0102473 0.0002784 36.80 <2e-16 ***
## test_set$FCH4 0.2612493 0.0109556 23.85 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01022 on 1800 degrees of freedom
## Multiple R-squared: 0.2401, Adjusted R-squared: 0.2396
## F-statistic: 568.6 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.089230 -0.004130 -0.000944 0.003110 0.124629
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063178 0.0001041 60.69 <2e-16 ***
## result$FCH4 0.5444722 0.0038757 140.48 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007653 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7326, Adjusted R-squared: 0.7325
## F-statistic: 1.974e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.084042 -0.006504 -0.001639 0.005549 0.038488
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0103947 0.0002558 40.63 <2e-16 ***
## test_set$FCH4 0.1863252 0.0082707 22.53 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009624 on 1800 degrees of freedom
## Multiple R-squared: 0.2199, Adjusted R-squared: 0.2195
## F-statistic: 507.5 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.202545 -0.004336 -0.000982 0.003647 0.075656
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0066319 0.0001111 59.68 <2e-16 ***
## result$FCH4 0.4825550 0.0041375 116.63 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00817 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.6538, Adjusted R-squared: 0.6537
## F-statistic: 1.36e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.054633 -0.006581 -0.001505 0.004878 0.062886
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0101317 0.0002792 36.28 <2e-16 ***
## test_set$FCH4 0.2488734 0.0103963 23.94 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01011 on 1800 degrees of freedom
## Multiple R-squared: 0.2415, Adjusted R-squared: 0.2411
## F-statistic: 573.1 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.142137 -0.003974 -0.000927 0.003076 0.110841
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0060593 0.0001046 57.94 <2e-16 ***
## result$FCH4 0.5414592 0.0038935 139.07 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007688 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7286, Adjusted R-squared: 0.7286
## F-statistic: 1.934e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.060266 -0.006622 -0.001702 0.005191 0.053183
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0105263 0.0002763 38.09 <2e-16 ***
## test_set$FCH4 0.2795353 0.0104466 26.76 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.0101 on 1800 degrees of freedom
## Multiple R-squared: 0.2846, Adjusted R-squared: 0.2842
## F-statistic: 716 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.138062 -0.004134 -0.001171 0.003135 0.087972
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063390 0.0001023 61.95 <2e-16 ***
## result$FCH4 0.5414730 0.0038100 142.12 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007523 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7371, Adjusted R-squared: 0.7371
## F-statistic: 2.02e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.108971 -0.006238 -0.001503 0.004592 0.062790
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0104501 0.0002724 38.37 <2e-16 ***
## test_set$FCH4 0.2390721 0.0101679 23.51 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01005 on 1800 degrees of freedom
## Multiple R-squared: 0.235, Adjusted R-squared: 0.2345
## F-statistic: 552.8 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.226219 -0.004132 -0.001025 0.003123 0.123061
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063691 0.0001059 60.16 <2e-16 ***
## result$FCH4 0.5330332 0.0039418 135.22 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007783 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7174, Adjusted R-squared: 0.7173
## F-statistic: 1.829e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.094774 -0.006536 -0.001607 0.004900 0.041756
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0101341 0.0002719 37.27 <2e-16 ***
## test_set$FCH4 0.2004374 0.0093309 21.48 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01016 on 1800 degrees of freedom
## Multiple R-squared: 0.204, Adjusted R-squared: 0.2036
## F-statistic: 461.4 on 1 and 1800 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.218365 -0.004080 -0.001016 0.003417 0.082211
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063685 0.0001111 57.31 <2e-16 ***
## result$FCH4 0.5090022 0.0041372 123.03 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008169 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.6775, Adjusted R-squared: 0.6775
## F-statistic: 1.514e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.068610 -0.006543 -0.001962 0.004819 0.084404
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0089865 0.0003133 28.68 <2e-16 ***
## test_set$FCH4 0.3193849 0.0122425 26.09 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01129 on 1800 degrees of freedom
## Multiple R-squared: 0.2744, Adjusted R-squared: 0.274
## F-statistic: 680.6 on 1 and 1800 DF, p-value: < 2.2e-16




## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.096717 -0.003831 -0.000957 0.003028 0.098357
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0054318 0.0001081 50.27 <2e-16 ***
## result$FCH4 0.5919832 0.0040230 147.15 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007944 on 7204 degrees of freedom
## (8631 observations deleted due to missingness)
## Multiple R-squared: 0.7504, Adjusted R-squared: 0.7503
## F-statistic: 2.165e+04 on 1 and 7204 DF, p-value: < 2.2e-16
## Warning: `as.tibble()` was deprecated in tibble 2.0.0.
## ℹ Please use `as_tibble()` instead.
## ℹ The signature and semantics have changed, see `?as_tibble`.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## StartDate EndDate AvailNum EffNum ObsNum AnnNEE_gCO2
## 1 2022-06-01 15:00:00 2023-04-27 13:00:00 15837 1477.062 9930 29.82144
## UncNEE_Rand_gCO2 UncNEE_uStar_gCO2 UncNEE_Comb_gCO2 AnnNEE_gC UncNEE_Rand_gC
## 1 63.53007 32.61282 71.41194 8.138726 17.33832
## UncNEE_uStar_gC UncNEE_Comb_gC
## 1 8.900535 19.48941
## StartDate EndDate AvailNum EffNum ObsNum AnnFCH4_gCH4
## 1 2022-06-01 15:00:00 2023-04-27 13:00:00 15837 2238.738 7034 6.101
## UncFCH4_Rand_gCH4 UncFCH4_uStar_gCH4 UncFCH4_Comb_gCH4 AnnFCH4_gC
## 1 0.2003837 0.1353286 0.2418004 4.568523
## UncFCH4_Rand_gC UncFCH4_uStar_gC UncFCH4_Comb_gC
## 1 0.1500504 0.1013362 0.1810639